Quantum phase transitions and dimensional reduction 
in antiferromagnets with inter-layer frustration 
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For magnets with a fuUy frustrated inter-layer interaction, we argue that the quantum phase 
transitions from a paramagnetic to an antiferromagnetic ground state, driven by pressure or magnetic 
field, are asymptotically three-dimensional, due to interaction-generated non-frustrated inter-layer 
couplings. However, the relevant crossover scale is tiny, such that two-dimensional behavior occurs in 
an experimentally relevant low-temperature regime. In the pressure-driven case the phase transition 
may split, in which case an Ising symmetry related to inter-layer bond order is broken before 
magnetism occurs. We discuss the relation of our results to recent experiments on BaCuSi2 06. 



In modern condensed matter physics, reduced dimen- 
sionality presents a fascinating avenue to novel effects 
arising from strong fluctuations. The standard realiza- 
tion of reduced dimensionality is via systems of chains 
or planes with a weak three-dimensional (3d) coupling in 
suitably structured materials. In such a situation, the 
physics is one-dimensional (Id) or two-dimensional (2d) 
at elevated energies (or temperatures), while it becomes 
3d in the low-energy limit. 

Recent experiments near magnetic quantum critical 
points (QCP) - both in heavy-fermion metals^"^ and in- 
sulating dimer magnets'* - have raised speculations about 
a rather different route to reduced dimensionality, namely 
through geometric frustration. The idea is that fully 
frustrated 3d interactions (being not necessarily weak) 
between 2d units effectively vanish in a well-defined low- 
energy limit, which in particular can be realized at a 
QCP. Then, it has been suggested that true 2d behavior 
could be observed at lowest energies."* 

Particularly interesting are the results on Mott- 
insulating quantum paramagnets, where the spin gap can 
be closed either by application of a magnetic field or 
pressure. The field-driven quantum phase transition 
(QPT), from a paramagnet to an XY-ordered antiferro- 
magnet (AF) at a field H — Hd, belongs to the universal- 
ity class of the dilute Bose gas, and the finite-temperature 
transition can be understood as Bose-Einstein condensa- 
tion of magnons. In BaCuSi206, the transition line has 
been found'* to follow Tc oc {H — Hd)'^ , with a shift 
exponent -0 = 1 characteristic of a 2d QCP, and indica- 
tions for a crossover to the 3d value ip — 2/3 at higher 
T. As BaCuSi2 06 has a body-centered tetragonal (bet) 
structure of Cu dimers, with a frustrated inter-layer cou- 
pling, the results have been interpreted as dimensional 
reduction at a QCP arising from geometric frustration. 

On the theory side, an investigation* of the ordered 
state of a frustrated double-layer Heisenberg model using 
spin-wave theory concluded that inter-layer order is sta- 
bilized through an order-from-disorder mechanism, thus 
rendering dimensional reduction ineffective. In fact, very 
similar results were obtained much earlier in spin-wave 
studies of frustrated bet magnets.^'*" However, the situa- 
tion is different at a QCP,''"** as the fluctuation-generated 



inter-layer coupling can be expected to be proportional 
to the square of the order parameter and hence absent 
at the QCP, supporting the interpretation of dimensional 
reduction from frustration. However, as we show below, 
higher-order processes turn out to be relevant in deter- 
mining the nature of the phase transitions. 

The purpose of this paper is an investigation of QPT 
in such Mott-insulating magnets with inter-layer frustra- 
tion, using a detailed symmetry analysis, field theory and 
bond-operator approaches. Our main results can be sum- 
marized as follows: (i) In general, geometric frustration 
cannot lead to dimensional reduction at asymptotically 
low energies, due to fluctuation-generated unfrustrated 
interactions. But their energy scale is tiny at quantum 
criticality, rendering 2d behavior observable even at very 
low temperatures, (ii) For the bet lattice, the ordered 
state breaks both the magnetic and an Ising symmetry. 
In the pressure-driven case, two scenarios are possible. 
Either there is a single transition breaking both sym- 
metries, or the transitions is split: the Ising symmetry 
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FIG. 1: (Color online) a) bet lattice of dimers. b) top view: 
two planes with sites shown as circles/crosses, and in-plane 
(inter-plane) coupling as solid (dashed) lines, c) Illustra- 
tion of the Z2 symmetry breaking, where the diagonal inter- 
plane bonds develop a spontaneous asymmetry (see text), 
together with the simplest lattice compatible distortion, d) 
T = phase diagrams for the coupled-dimer model (1), for 
the pressure-driven case with and without splitting (see text) 
and the field-driven case. 
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is broken first, which relieves the intcr-laycr frustration, 
and the subsequent magnetic transition is conventional, 
Fig. 1. 

Model Hamiltonian. To be specific, we will consider 
a coupled-dimer system on a bet lattice (but we believe 
our results to be more general), with the Hamiltonian 

W — J Sin\ • Sin2 + J ^ ^ Sinm'Sjnm (1) 

in {ij)nm 

I \ ^ jmm' c q fj \ ^ q 

iAnrnm' inm 

where m = 1,2 labels the two spins 1/2 of each dimer, 
i,j are the dimer site indices in each layer, and n is 
the layer index. J and J' arc the AF intra-dimer and 
in-plane inter-dimer couplings, respectively, while J™™ 
represent the frustrated inter-layer couplings (with some 
specific inter-dimer structure given by the mm' depen- 
dence). The runs over four sites such that the sites 
{in) and (t + A, n + 1) arc nearest neighbors in z direc- 
tion. This Hamiltonian is assumed to be relevant for the 
material BaCuSi206^ (neglecting here the orthorhombic 
distortions in the low-temperature phase^^). 

Phases. For J ^ J', | J^], the zero- field ground state of 
W is a paramagnetic singlet, with gapped triplet excita- 
tions. If J' dominates, an AF phase with in-plane order- 
ing wavevector Q = (tt, tt) is established, and the order 
in z direction is frustrated, as discussed below. (For large 
\Jz \ the in-plane order is ferromagnetic.) Applying a field 
to the large- J quantum paramagnet leads to a Zeeman 
splitting of the triplet excitations. At a critical field Hd, 
the gap of the lowest mode closes, and a QPT to a gaplcss 
canted phase occurs. Upon further increasing the field, 
the system is driven into a fully polarized state at i?c2- 
The phase diagram at T = is thus similar to that of 
the much-studied bilayer Heisenbcrg model. ^^'^^ 

AF phase: Order from disorder. While AF ordered 
planes of classical moments on the bet lattice are uncou- 
pled, zero-point fluctuations of quantum spins lift this 
large degeneracy. A spin-wave calciilation, for spins S 
with couplings J', Jz and a helical order with wavevec- 
tor (7r,7r, Qz), yields an inter-layer contribution to the 
groimd-state energy of the form {—J^S/,J'cos^Qz), fa- 
voring collinear (i.e. parallel or antiparallel) order be- 
tween adjacent planes. ^'^ In addition, higher-order terms, 
in a calculation for a more general ordering pattern, ac- 
tually stabilize AF order between .2nc?-neighbor planes. 
Thus, fully 3d order is stabilized within the stacks of 
"even" and "odd" planes, but a residual Z2 degeneracy 
is left intact, corresponding to a spin inversion in every 
second plane this represents a true symmetry of the AF 
on the bet lattice, and will be spontaneously broken in 
the ordered phases. 

Symmetries and magnetic order parameter. The in- 
plane magnetism is unfrustrated, allowing to define an 
order parameter (^„(r||), where f|| is the in-plane coordi- 
nate and n the layer index, with the local magnetization 
operator given by m„(f^|) = exp(i(5 • r||)<?!)n(r||). Expand- 



ing the tight-binding dispersion on the bet lattice aroimd 
Q, i.e. taking the in-plane continuum limit, we arrive at 
a (f)^ theory for the magnetic fluctuations: 

= I dTd^k\\^Y'^^+^^^)^n{h) 



<^dyn 



(2) 



where ^i>n(A;||) is the real order-parameter field, fcy = 
now corresponds to physical in-plane momentum (7r,7r), 
c is a velocity, and rj represents the spatial anisotropy 
(i.e. ?7 ~ Jz/J'). Further, iS^dyn encodes the dynamics 
of the spin fluctuations (which depends on whether an 
external field is present or absent), and 5^4 is the local 

quartic self-interaction, (5^4 = uq J drd^r|| X]„[0^(?^|)]^- 
Effects of the geometric frustration completely dominate 
the inter-layer dispersion: it is not only suppressed at 
zero wavevector, but also acquires a sign alternating from 
quadrant to quadrant away from it.^ 

The action (2) represents the "bare" order-parameter 
dynamics; additional symmetry-allowed terms will be 
generated upon integrating out high-energy modes. We 
analyze the symmetries of the problem: In the zero-field 
paramagnetic phase, time and space inversion as well as 
SU(2) spin rotations are unbroken. The system is also 
invariant under 90-degree in-plane rotations, however, 
the frustrated geometry dictates that this is accompa- 
nied by a relative sign change of the order parameter in 
two neighboring planes: 



ky , ky > kx , (j^n ^ ( l)"^n • 



(3) 



Obviously, the action S0 (2) is invariant under this trans- 
formation. Now we can list the lowest-order additional 
terms allowed and forbidden by symmetry: 

Allowed: ??'(?„-0„+2, 'iti(0„-0„+i)^, M20^(^+i, •■(4) 
Forbidden: ^„ • $n+i, kx$n ■ $n+i,-- (5) 

Thus, single-particle hopping between neighboring planes 
must be suppressed at least by a factor k^ky (i.e. it is 
frustrated), as in (2). However, single-particle hopping 
between 2nd-neighbor planes is allowed even at fcy = 0, 
as are density interactions between and (pn+i- 

The first two terms in (4) are of crucial importance 
for the following analysis; they obviously correspond to 
additional terms to the microscopic Hamiltonian (1) 

(6) 

in'mm' 

Wcoll — ^ ^ "^coll i^inm ' Si-\-A,n-\-l,m' r, (7) 

iAnrnm' 

which can be present in any real material. However, even 
if absent in the "bare" model, they will be perturba- 

tivcly generated. The corresponding lowest-order dia- 
grams, quadratic in the interaction vertex uq, are shown 
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FIG. 2: (Color online) Diagrams relevant for the order- 
parameter theory S^. Lines are (j)„{k^^) propagators, the circle 
is the local four-point vertex (oc uo), the cross is the inter- 
layer hopping (cx rj). a) Unfrustrated vertical 2nd- neighbor 
coupling - this is the ry' term in (4) which is responsible for 
3d behavior at lowest energies, b) Inter-layer density inter- 
action, generating the ui and U2 terms in (4). ui < leads 
to collinear inter-plane spin correlations, c) Vertical near- 
neighbor coupling in the ordered phase (the open square de- 
notes the condensate). 

in Fig. 2a, b. Even if the u terms are (marginally) irrele- 
vant in the RG sense at the QCP of the theory (2) (which 
is the case at the Hd critical point), the 77'</)„ • (f>n+2 term 
is strongly relevant. 

Thus we arrive at a central result: There exists an 
energy scale in the 4> dispersion, Ez ^ rj' , corresponding 
to an unfrustrated 3d coupling in the paramagnetic phase, 
which is finite also at the QCP. Then, the behavior of 
observables at energies or temperatures below will 
be fully 3d (with a dimensional crossover to 2d behavior 
above i?z)j and no true dimensional reduction occurs. 

A straightforward calculation of the diagram in Fig. 2b 
gives a negative coefficient ui oc —rf' for the {(pn ■ (pn+if 
interaction, thus stabilizing collinear spin correlations be- 
tween adjacent planes, consistent with the spin-wave re- 
sult for the ordered phase. In fact, an effective Hamil- 
tonian of the form (7), with Jcoii oc — J|/(J'S''^), was 
proposed before in* to mimic the collinear order-from- 
disorder mechanism within linear spin-wave theory. 

Ising order parameter. As emphasized above, the mag- 
netically ordered phases break an additional Z2 symme- 
try - in the (j) language this degree of freedom corresponds 
to the sign of (0„) • (0„+i). It is useful to introduce an as- 
sociated (inter-plane) Ising order parameter ^n+i/2{r\\)- 
It transforms as a singlet under spin rotations, hence 
symmetry dictates that can be described by a stan- 
dard theory iSip, similar to (2), but with mass to^ 
and an unfrustrated Yertical hopping 'I'„-i/2^n+i/2- The 
physical content of 5* becomes clear from its interaction 
with (f), which is of Yukawa type: 

^0* = A / drd^kii 'I'n+1/2 $n ■ (fn+i (8) 

n 

with A a coupling constant. (Additional couplings '^(f'n 
do not modify the physics to be discussed below.) 

One vs. two transitions. The full theory iS^-l-iS^-l-iS^* 
admits two distinct scenarios: (A) A single transition 
driven by the condensation of (j) - here, the coupling A 
generates a non-zero expectation value for ^' as well, be- 
cause the ui term in leads to non-zero (0„ • 0n+i) 
(Fig. 2c). (B) Two transitions: First, "if condenses, which 
only modifies the quadratic part of the 4> action, and (j) 



orders in a second, subsequent transition. 

Scenario (B) is appealing: From (8) it can be seen 
that a ^' condensate generates an unfrustrated vertical 
hopping of 4> through the term A(5')(/)„ • 4>n+i- Thus, the 

ordering transition removes the magnetic inter-plane 
frustration. In the lattice model, this can be understood 
as spontaneous bond order which modulates the vertical 
magnetic couplings within each unit cell and can easily 
couple to lattice distortions. Fig. Ic, i.e., ^E" ordering is a 
structural phase transition. 

In zero field, a model with negative ui may follow sce- 
nario (B) for the pressure-driven transition: 5* can be un- 
derstood as composed of two 4> quanta, hence ~ 2m0. 
The Ml term mediates an attraction between and 
4'n+i- Approaching the (j) ordering transition, m>ii can 
become smaller than m^, implying that 5" condenses be- 
fore (f>. This requires a sufficiently strong (otherwise 
no true two-particle bound state is generated, and the 
transition remains in scenario (A)). 

Lattice theory. We have studied the coupled-dimer 
model (1) using the bond-operator approach. Starting 
from a singlet product state on dimcr sites i, Yii Njs), 
we define bosonic operators tl^ which create local triplet 
excitations \i. a) — tj^ \i. s) where a = 0, — and 
\h+) = -ITT), K,-> = Ui) and |z,0) = (iTi) + 
|iT>)/\/2- The Hamiltonian (1) can be re-written in 
triplet operators. "'^'^"^^ The quadratic part reads: 

13 ^ 

n2 = J2 - 4" ^9" + -f (Vai-g-a + h.c.) (9) 

qa 

with a — —a, — J + B^, = J'(cosg^ -I- cosqy) + 
2J2 cos((7^/2) cos(gy/2) cosg^, where = J]^ + J^^ — 
J^^ — J^^, and the field H is in z direction. In addition 
to Ti.2 the Hamiltonian contains a quartic in-plane triplet 
term ^411;" for .h^ = J^" + + J]^ + 7^ and 
J3± = Jii - J22 ^ (•ji2 _ j2i) ^ Q quartic and cubic 
inter-plane interactions, Ti^z and Ti^z, arise, respectively. 

The eigenvalues of 7^2, (^qa = - Bl)^/"^ - aH, are 
independent of Qz at in-plane wavevector — {tt,tt). 
In this linearized bond-operator theory, interactions be- 
tween the order-parameter fluctuations - represented by 
triplon quasiparticles - are ignored. The most impor- 
tant corrections arise from the hard-core constraint Tiu, 
J2a^la^ia — 1-^^ FoUowing the self-consistent diagram- 
matic approach of Ref. 14, we have calculated the renor- 
malized triplon dispersion, focusing on how interactions 
lift the degenerate vertical spectrum, Fig. 3. Taking into 
account Hu yields a dispersion along (tt, tt, g^) propor- 
tional to JzCos{2qz), corresponding to the process in 
Fig. 2a. Further analysis shows that a similar term, but 
with prefactor J^, is generated to second order in 7^42. 
(This cannot be read off from the order-parameter the- 
ory because of the assumed local form of the quartic 
term.) Hence, the bandwidth along (7r,7r,gz) is at most 
of order J^.^* (This result is not changed by Tisz.) 

To assess the dynamics of the Ising parameter ^, we 
have studied bound states of two triplons in the singlet 
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FIG. 3: (Color online) Triplon dispersion (solid) and bound- 
ary of the two-particle continuum (shaded), calculated using 
bond operators, for J' /J = 0.15, = = 0.2J', H = 0. 
The inset show an energy zoom into the almost flat vertical 
dispersion. The dominant effect of is a Zeeman shift of 
—aH of the triplons. 



channel by solving a Bethe-Salpeter equation, with a 
four-point vertex as input. To avoid the necessity for full 
self-consistency in the two-particle sector, we have added 
to Hi by hand the biquadratic TicoU (7). The effective 
binding force is then given by (J^^ — Jco\i)] note that the 
binding from 7^42 is due to a subtle quantum effect of 
singlet formation, which is not contained in the order- 
parameter field theory, but is in fact very common for 
frustrated spin-1/2 systems. "'^^'^^ 

For sufficiently strong positive (J42 — Jcou) we find 
a bound state below the two-particle continuum, whose 
wavefunction changes sign under 90-degree in-plane rota- 
tions of the internal coordinate, consistent with (3) and 
Fig. Ic. Close to the pressure-driven phase transition the 
bound state (with dispersion minimum at total momen- 
tum Q = 0) is pulled below the sm(7fe-particle gap, due 
to the weak vertical triplon dispersion. Clearly, a con- 
densation of this bound state corresponds to the Z2 sym- 
metry breaking advocated above. This existence of the 
bound state depends on microscopic details, and a com- 
prehensive numerical analysis is difficult due to finite-size 



effects. (For the parameter values of Fig. 3, our calcula- 
tions indicate no bound state.) 

Pressure vs. field tuning. So far, most considerations 
were for H — 0. In finite field, the spin symmetry is re- 
duced to U(l), and the expressions in (2,4,5,8) are mod- 
ified accordingly. However, the symmetry analysis in the 
paramagnetic phase, leading to (4,5), remains valid. The 
most important difference is in the bound-state behav- 
ior: Close to Hci a possible singlet bound state involves 
at least one high-energy triplet and will never condense. 
Then, scenario (A) always applies, with a 3d QPT which 
breaks the U(l) x Z2 symmetry and obeys mean-field 
exponents. 

Discussion. We conclude that dimensional reduction 
in quantum-critical frustrated bet magnets does not oc- 
cur at lowest energies, due to interaction effects. How- 
ever, the crossover scale below which 3d behavior is es- 
tablished (which is given by the vertical dispersion) is 
tiny, Ez oc / (instead of Jz as in an unfrustrated 
system). Above E^, the shift exponent ip will take its 2d 
value. In the pressure-driven case, depending on micro- 
scopic parameters the ordering transition may be split, 
then bond order occurs before magnetic order. 

Relevance to BaCuSi2 0Q. Recent neutron scattering^*^ 
hints on the presence of multiple triplon excitations in 
zero field, which indicate an enlarged unit cell with in- 
equivalent dimers. Assuming that the lattice modulation 
occurs along the c axis, one arrives at a qualitatively 
distinct scenario for "dimensional reduction" : The con- 
densate established at Hd is strongly inhomogeneous in z 
direction, and hence effectively 2d - note that frustration 
is not a required ingredient in this scenario. (This conclu- 
sion is supported by the fact that the measured vertical 
triplon dispersions are tiny for all in-plane wavevectors,^^ 
at variance with the results in Fig. 3.) 
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